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We analyze two-center interference effects in the yields of ionization of a dissociating hydrogen 
molecular ion by an ultrashort VUV laser pulse. To this end, we performed numerical simulations of 
the time- dependent Schrodinger equation for a H J model ion interacting with two time-delayed laser 
pulses. The scenario considered corresponds to a pump-probe scheme, in which the first (pump) 
pulse excites the molecular ion to the first excited dissociative state and the second (probe) pulse 
ionizes the electron as the ion dissociates. The results of our numerical simulations for the ionization 
yield as a function of the time delay between the two pulses exhibit characteristic oscillations due to 
interferences between the partial electron waves emerging from the two protons in the dissociating 
hydrogen molecular ion. We show that the photon energy of the pump pulse should be in resonance 
with the ag — Gu transition and the pump pulse duration should not exceed 5 fs in order to generate 
a well confined nuclear wavepacket. The spreading of the nuclear wavepacket during the dissociation 
is found to cause a decrease of the amplitudes of the oscillations as the time delay increases. We 
develop an analytical model to fit the oscillations and show how dynamic information about the 
nuclear wavepacket, namely velocity, mean internuclear distance and spreading, can be retrieved 
from the oscillations. The predictions of the analytical model are tested well against the results of 
our numerical simulations. 



PACS numbers: 33.80.Rv, 42. 25. Hz, 42. 50. Hz 

I. INTRODUCTION 

Double-slit like interferences similar to those observed 
by Young in his experiment with light [1 occur also in 
the (photo-) ionization of diatomic molecules. Concep- 
tually, the partial electron waves ejected from the two 
atomic centers of the molecule take the role of the coher- 
ent light waves emerging from the two holes in Young's 
experiment. The phenomenon was first discussed by Co- 
hen and Fano to explain the observation of oscilla- 
tions in the total photionization cross section of the H2 
molecule as a function of the photon energy [3^ . Accord- 
ing to their analysis, the interference pattern depends 
periodically on the ratio of the internuclear distance to 
the wavelength of the photoelectron and is superimposed 
on the quickly decreasing photoabsorption cross section. 
In recent years. Young- type interferences have been ob- 
served experimentally and studied theoretically for the 
photoemission of (core) electrons from H2" and H2 |4HTQ] 
and heavier molecules pJT[T4] . Furthermore, interference 
effects in the ionization of molecules by heavy ion [15lfT6] 
and electron impact [171 [IH] have been investigated. Re- 
cently, the idea of Cohen and Fano was extended to the 
double photoionization of H2 [TOHSS], the photoioniza- 
tion of more complex molecules [24], and laser induced 
multiphoton ionization of molecules [25] [26] . 

In photoionization of a diatomic molecule the double- 
slit interferences primarily depend on the momentum of 
the photoelectron, /c, the internuclear distance between 
the two atoms, i?, and the angle between the emission 
direction of the electron and the molecular axis. In 
past work, photoionization of molecules from the ground 
electronic state and corresponding interference patterns 
in the angular distributions of the photoelectron in the 



molecular body fixed frame was studied. Even the small 
variations of the internuclear distance in the vibrational 
states of a bound molecule could be observed in the pat- 
terns [mils]. 

However, the dependence of the interference patterns 
on the internuclear distance becomes more apparent in 
the photoionization from dissociating molecules. In a 
recent set of experiments negative iodine molecules I^ 
were excited by a pump laser pulse to a dissociative 
state, in which the molecule fragments into a negative 
ion I~ and a neutral atom I [27H29). During the disso- 
ciation the molecule was ionized by a second pump laser 
pulse and the photoelectron spectra were measured. The 
anisotropy parameter obtained from the photoelectron 
angular distributions was found to oscillate as a func- 
tion of the pump-probe time delay, in other words as 
a function of the increasing internuclear distance in the 
dissociating molecule. As in the previous examples, the 
oscillation arises due to the interference between the par- 
tial electron waves emerging from the two fragments. 

Very recently, the probe step in such kind of pump- 
probe experiments, namely the photoionization of the 
dissociating molecule, has been studied theoretically us- 
ing first-order perturbation theory for the prototype ex- 
ample of a II2" molecular ion, where the positions of the 
nuclei were fixed at different internuclear separations [30] . 
Assuming a nuclear wave packet of Gaussian shape with 
a certain but fixed width Ai? centered at an internu- 
clear distance Rq^ the angular and momentum distribu- 
tion of the photoelectron emitted from the first excited 
electronic state (cr^^-orbital) of II2" was calculated. The 
two-slit interferences were explored for different orienta- 
tions of the molecular axis with respect to the polariza- 
tion direction of the probe laser. 
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FIG. 1: (Color online) Scheme of the pump-probe scenario, 
in which a first (pump) pulse initiates the dissociation by 
excitation of a wavepacket to the au-state of H^. With a 
time delay r a second (probe) pulse ionizes the dissociating 
molecular ion. The green curve denotes the energy of two 
protons as a function of internuclear distance. 



In the present work we complement this recent theo- 
retical study by results of numerical simulations in which 
we investigate the complete pump-probe scenario accord- 
ing to the following scheme (see Fig. 1 for illustration). 
We consider a H2" molecular ion initially prepared in its 
electronic (cr^) and vibrational {u = 0) ground state. The 
interaction with an ultrashort pump pulse resonantly ex- 
cites part of the wave packet into the repulsive a^-state, 
thus initiating the dissociation of the molecule. The 
ungerade state corresponds to an antisymmetric elec- 
tronic wave function consisting of a coherent superpo- 
sition of atomic orbitals localized at the two dissociating 
protons, which establishes a double-slit scenario. During 
the dissociation the protons first accelerate in the steep 
part of the molecular ion potential and then propagate 
with almost constant velocity. Furthermore, the width 
of the nuclear wave packet increases as a function of in- 
ternuclear distance due to dispersion. This requires to 
generate a well-confined nuclear wave packet in the exci- 
tation (pump) step in order to resolve the interferences 
in the ionization yields, as we will show below. In our 
simulations, we then modeled the interaction of the dis- 
sociating molecular ion with a second ultrashort probe 
laser pulse to determine the ionization probability as a 
function of the time delay between the dissociating and 
ionizing pulses. As we will show below, our results ex- 
hibit the expected Young- type interferences. We finally 
develop an analytical model in order to use such obser- 
vations as an imaging tool for the velocity, internuclear 
distance and the spread of the nuclear wave packet. 

The paper is organized as follows. In section [n| we in- 
troduce the H2" model ion used to perform the numerical 
simulations. Then, in section III we study the dissoci- 



cuss the analogy of the probability density distribution 
of the dissociating system to a two-slit scenario for differ- 
ent parameters of the pump pulse. Next, results for the 
ionization yields will be presented in section |IV[ where 
the double-slit interference patterns will be identified and 
analyzed in view of imaging the characteristic features of 
the nuclear wave packet. Finally, we conclude with a 
short summary. 



II. H+ MODEL FOR NUMERICAL 
SIMULATIONS 

For our analysis of the ionization of the dissociating hy- 
drogen molecular ion following the pump-probe scheme, 
illustrated in Fig. 1, we make use of a model for H2" 
in which the electronic and nuclear motion are restricted 
to one dimension. The model accounts for the interac- 
tion of the electron with the external fields as well as for 
non-Born-Oppenheimer effects. In the model, the polar- 
ization direction of the two linearly polarized pump and 
probe laser pulses coincide with the internuclear axis of 
the molecular ion as well as with the restricted electronic 
motion. The field-free Hamiltonian of this model system 
is given by (we use Hartree atomic units, e = m = h = 1) 
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where P = —id/dR and i?, and p = —id/dz and z are 
the linear momentum operators and the positions of the 
relative nuclear (protons) and electronic coordinate, re- 
spectively. fiR = M/2 and /i^ = 2M/(M + 1) are the 
reduced masses where M = 1836 a.u. is the mass of the 
proton. We have chosen soft-core parameters = 1 and 
ap = 0.03 [32]. The interaction of the electrons is taken 
into account in dipole approximation and velocity gauge 
by: 
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where Aq^i is amplitude of the vector potential, uji is the 
frequency of the carrier wave, Ni and Ti are the number 
of cycles and the period of the carrier wave of either the 
dissociating (/ = d) or the ionizing (/ = i) pulse, respec- 
tively, and r is the time delay between the two pulses. 

The corresponding time-dependent Schrodinger equa- 
tion 



ation of the H2" model ion by the pump pulse and dis- 
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was solved using the Crank-Nicolson method [33^. The 
sizes of the simulation box were 1 to 37 a.u. and -409.6 
a.u. to 409.6 a.u. in R and z direction, respectively, 
with a grid spacing of Ai?=0.03 a.u. and Az=0.1 a.u.. 
The time step in the simulations was At=0.02 a.u.. The 
electronic-vibrational ground state of the model hydro- 
gen molecular ion was obtained using imaginary time 
propagation, the ground state energy was -0.77 a.u., the 
equilibirum distance 2.6 a.u. and the ionization potential 
for ^ oo about 0.71 a.u.. 

We have chosen to restrict the electronic motion to one 
dimension along the polarization axis and the internu- 
clear axis in our model, since we analyzed the ionization 
of the dissociating hydrogen molecular ion up to rather 
large internuclear distances (37 a.u.) and used a largely 
extended grid in z-direction (> 800 a.u.) in order to keep 
the wavefunction on the gird for further analysis. In re- 
cent studies, higher-dimensional calculations have been 
performed, in which the full three-dimensional electronic 
motion [10] or the two-dimensional electronic motion in 
cylindrical coordinates have been taken into account [34]- 
[37] . These higher-dimensional calculations are usually 
either performed on smaller grids or make use of high 
performance supercomputers. Here, we are interested in 
the analysis of the double-slit interference effects in the 
ionization signal from the dissociating molecular ion as 
a function of the internuclear distance R. These effects 
are expected to be present independent of the emission 
direction of the electron, which justifies a restriction of 
the electronic motion to one dimension. 



III. DISSOCIATION DYNAMICS 

First, we analyze the dissociation dynamics of the hy- 
drogen molecular ion after the interaction with the pump 
laser pulse. To this end, we performed numerical simu- 
lations of the model hydrogen molecular ion interacting 
with just one ultrashort laser pulse. We filtered out the 
electronic-vibrational ground state after the interaction 
and monitored the wavepacket in the excited state as a 
function of time. In Fig. [2] we compare two snapshots 
of the temporal evolution of the excited state probability 
density as a function of z and R (panels (b) and (c)) with 
the corresponding distribution in the ground state (panel 
(a)). The calculations were performed for a dissociating 
pump pulse with ujd = 0.38 a.u. (period = 16.22 a.u.), 
corresponding to a central wavelength of 117 nm, which 
is resonant with the cig — du transition in our model. The 
full pulse contained three cycles {Ni = 3), i.e. a full pulse 
length of 48.66 a.u., about 1.2 fs. The amplitude of the 
vector potential was Aq^^ = 0.04 a.u., corresponding to a 
peak intensity of 8.4 x 10^^ W/cm^. At the end of the 
pulse, the total excitation probability, i.e. the probability 
after projecting out the ag ground state, was 2.26%. 

Fig. |2]shows, as expected, that the excited state proba- 
bility density propagates to larger internuclear distances 
after the interaction with the pump pulse, which indi- 
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FIG. 2: (Color online) Comparison of the electron probabil- 
ity density in the ground state of the model (panel (a)) 
with the temporal evolution of the density in the first excited 
state of hydrogen molecular ion following an excitation of the 
molecular ion by an ultrashort laser pulse with a central wave- 
length of 117 nm and a full pulse length of 48.66 a.u. (about 
1.2 fs). Shown are snapshots of the excited state wavepacket 
density at (b) 30Td (11.782 fs) and (c) GOTd (23.564 fs). 
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FIG. 3: (Color online) (a) Expectation value of the internu- 
clear distance as a function of time, (b) Enlarged view of (a) 
for a time interval of 2 fs (16 fs to 18 fs). For the sake of com- 
parison, we have also plotted the electric field of a 10-cycle 
pulse at a wavelength of 45 nm, which was used as a probe 
pulse below. 
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FIG. 4: (Color online) Electron probability density (after re- 
moval of the electronic- vibrational ground state) following the 
interaction with pump laser pulses at 117 nm having a full 
pulse length of (a) 5 fs at t = 61Td (23.912 fs) and (b) 20 fs 
at t = 80Td (31.360 fs). 



cates a dissociation of the molecular ion. One sees that, 
at a given internuclear distance, the density is almost 
equally (de-) located at both protons and, hence, it forms 
the anticipated set-up, similar to a double slit. We note 
that the initially very well confined nuclear wave packet 
with a narrow distribution in R does spread significantly 
due to dispersion as the time proceeds. 

Next, we present in Fig. [3] the result from the nu- 
merical simulations for the expectation value of the in- 
ternuclear distance {R) as a function of time after the 
interaction with the pump pulse. Please note that the ex- 
pectation value was calculated after the remaining popu- 
lation in the ground state was removed. It is clearly seen 
that the wave packet gets accelerated over the first few 
femtoseconds due to the steep gradient of the repulsive 
potential energy curve for the cr^^-state at short internu- 
clear distances. However, we observe that about 4-5 fs 
after the end of the interaction with the pump pulse the 
velocity of the wave packet is almost constant and (R) 
increases linearly with time. From the expanded view 
in Fig. [sjb) we notice that the expectation value (R) 
changes by 0.86 a.u. (±0.02 a.u.) within 1 fs in the lin- 
ear regime, i.e. the nuclear wave packet propagates with 
a velocity of 0.0104 a.u. (±0.0003 a.u.). 

We also studied the dissociation dynamics for pump 
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FIG. 5: (Color online) Electron probability density (after re- 
moval of the electronic-vibrational ground state) following the 
interaction with pump laser pulses at 800 nm with a peak in- 
tensity around 10^"^ W/cm^ and having a full pulse length of 
(a) 1 cycle (i.e. 2.7 fs) at 23.564 fs and (b) 3 cycles (i.e. 8.1 
fs) at 25.872 fs. 



pulses with longer pulse lengths as well as with differ- 
ent frequencies in order explore the parameter regime of 
the pump pulse for creating a well-confined two-center 
scenario in the dissociating molecular ion. In Fig. [4] we 
present the excited state probability density generated 
by different pump pulses of 5 and 20 fs, which are reso- 
nant with the (jg — cFu transition. While the wavepacket 
generated with the 5 fs pulse is still well confined at large 
internuclear distance, we observe that the spread of the 
wavepacket is significantly larger after excitation with a 
20 fs pulse. This is due to the fact that for longer pulses, 
the (effective) interaction time between the molecular ion 
and the pump pulse increases. Consequently, the width 
of the initial wavepacket created in the excitation step 
increases, which leads to a wider wavepacket at larger in- 
ternuclear distances. Please note that both these pulses 
create a significant amount of vibrational excitation in 
the bound molecular ion (cf. density at small internu- 
clear distances in Fig. |4|. 

For the sake of comparison, we show in Fig. [5] the cor- 
responding excited state probability densities generated 
by the interaction with a non-resonant pump laser pulse 
operating at 800 nm and total pulse lengths of (a) 1 cy- 
cle (2.7 fs) and (b) 3 cycles (8.1 fs). At this wavelength 
it is well-known that the dissociation of H2" occurs via 
other mechanisms such as above threshold dissociation 
j38] , bond softening [39l [40] , or bond hardening [41] . In- 
stead of one (narrow or wide) wavepacket, as in the case 
of a resonant one-photon excitation with the VUV pulses, 
we observe a train of nuclear wave packets for both pulse 
lengths. Also, in this case the vibrational states are much 
more excited than in case of an excitation by the VUV 
pulses. 

We may summarize our results so far in view of our 
goal to observe double-slit interference effects from the 
dissociating hydrogen molecular ion. For this goal, we 
require a single well-confined nuclear wavepacket, which 
corresponds to a narrow distribution of the internuclear 
distances (which shows certain similarities with a double- 
slit having a well-defined separation of the slits). From 
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the results presented above, we conclude that the sys- 
tem should be prepared in the vibrational ground state, 
the pump pulse should be in resonance with the Gg — 
transition (in order to avoid trains of wavepackets as in 
the case of dissociation driven by the 800 nm pulses) 
and as short as possible (in order to confine the initial 
wavepacket as much as possible). Furthermore, we infer 
from our result for the velocity of the nuclear wavepacket 
(cf. Fig. |3jb)) that the ionizing probe pulse should be 
less than about 2 fs. Otherwise the wavepacket sweeps 
over a rather large region of internuclear distances dur- 
ing the interaction with the probe pulse which will likely 
smear out any interference patterns in the ionization sig- 
nal, as long as the internuclear distance at the moment 
of ionization is not detected in the experiment. 

IV. PHOTOIONIZATION OF DISSOCIATING 
H+ MOLECULE 

In this section we present the results of numerical 
simulations in which we took the complete pump-probe 
scheme into account (Fig. [T]). First, we focus on the 
identification of interference patterns in the yields of ion- 
ization of the dissociating hydrogen molecular ion by a 
time-delayed probe pulse and then we analyze the re- 
sults in view of an imaging of the nuclear wave packet. 
For our analysis we used an ultrashort dissociating pump 
pulse, resonant with the (Jg — cFu transition (co'c^=0.38 a.u., 
A=117 nm) and a full pulse length of 48.66 a.u. (1.2 fs, or 
3-cycles) and a peak intensity of 8.4 x 10^^ W/cm^. The 
corresponding dissociation dynamics of the model H2" is 
represented by the results shown in Figs. |2]and[3] In each 
simulation, we project out the ground state of the H2" 
model ion after the interaction with the pump pulse, as- 
suming that any ionization signal from the bound molec- 
ular ion can be well separated in the experiment. 

A. Double-slit interferences 

We performed a series of numerical simulations in 
which we considered the application of a second photo- 
ionizing probe pulse as a function of time delay r. After 
the end of the probe pulse we calculated the ionization 
yield via the integral over the outgoing parts of the elec- 
tron wave packet as (please note, that we projected out 
the ground state of the H2" model after the pump step): 

00 00 

P+=lim f dR I dz\^{z,R,t)\'^ . (5) 

t^OO J J 

|2| = 10+f 

We continued the propagation of the wavefunction un- 
til the ionization signal determined on the grid did not 
change as a function of time. 

According to previous analysis the photoionization 
yield from a diatomic molecule is expected to be mod- 
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FIG. 6: (Color online) Results of numerical simulations for 
the ionization yield using two different wavelengths for the 
ionizing pulse: (a) 22 nm and (b) 45 nm. 

ulated as [2l[T3]: 

P+(k, R) ocl + cos(k • R + ^(k, R)) , (6) 

with <l>(k, R) = (/)orb + 0scat(k, R), where 0orb accounts 
for the phase of the symmetry of the orbital from which 
the electron is emitted (0orb = 7r/2 for the present case, 
as the electron is in the a^^-orbital), and (/)scat(k, R) rep- 
resents a phase acquired due to the scattering of the 
electron by the neighboring proton [13 . Since the sym- 
metry of the electronic molecular wavefunction does not 
change during the dissociation, (/)orb does not depend on 
the internuclear distance and the photoelectron momen- 
tum. In contrast, one may expect that in general the 
acquired scattering phase depends on both these param- 
eters. Since in our calculations we restricted the electron 
motion to one dimension along the internuclear axis of 
the molecular ion, Eq. (|6| reduces to: 

P^{k,R) (xl^cos{kR^^{k,R)) . (7) 

We considered two different ionizing pulses. The first 
one had a frequency of uji = 2.07 a.u. (56.33 eV per 
photon), corresponding to a wavelength of 22 nm and 
a period of Ti = 3.03 a.u. (73 as), Ni = 10 (i.e., full 
pulse length of 730 as) and Ai = 0.03 a.u. (i.e., peak 
intensity of 1.35 x 10^"^ W/cm^). The parameters of the 
second pulse were uJi = 1.01 a.u. (27.54 eV per photon), 
corresponding to a wavelength of 45 nm and a period 
of Ti = 6.20 a.u. (150 as). The number of cycles was 
Ni = 10 (i.e., full pulse length of 1.5 fs) and the amplitude 
of the vector potential was Ai = 0.06 a.u. (i.e., peak 
intensity of 1.3 x 10^"^ W/cm^). 

In Fig. |6] we show the numerical results for the ioniza- 
tion yield as a function of the pump-probe time delay 
for the two ionizing pulses. We clearly observe oscil- 
lations of the ionization yields due to the interferences 
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of the partial electron wave packets. The amplitudes 
of the oscillations decrease, which we attribute to the 
spreading of the nuclear wavepacket (for further analy- 
sis, see section IV. B.). Over the range of time delays 
studied, the frequency of the oscillations remains prac- 
tically constant (0.020 ± 0.001 a.u. for ionization at 45 
nm, and 0.038 ± 0.005 a.u. at 22 nm). Since the ion- 
ization potential of the dissociating hydrogen molecule 
in the present range of internuclear distance does vary 
very little and, hence, the photoelectron momentum re- 
mains almost constant, we conclude that the phase ^ 
and, hence, the scattering phase (/)scat is independent of 
the internuclear distance R (or, of the time delay r). 
This result makes sense, since in the present model we 
restricted the electron motion to one dimension along the 
internuclear distance and at large internuclear distances 
the Coulomb potentials located at the two protons are 
well separated and vary just slightly with an increase of 
R. Thus, the scattering process of the electron at the 
neighboring Coulomb potential is (almost) independent 
of the value of R and, hence, ^{k^R) = ^{k) for large R 
in our model. 

Before we proceed, we may mention that we have also 
considered probing the dissociating molecular ion using 
near-infrared laser pulses at 800 nm. Due to the non- 
perturbative coupling between the electron and field at 
this wavelength, the ionization of the electron proceeds 
via tunneling or multiphoton above-threshold ionization. 
Consequently, the electron is not emitted at a given en- 
ergy (within a certain width) but with a range of electron 
energies that is rather broad. In order to observe any 
trace of the double-slit interference as a function of the 
internuclear distance, it is therefore necessary to restrict 
the electron energy in an experiment as well as in the 
numerical simulations. In other words, instead of total 
ionization yields, we have to obtain partial differential 
ionization yields 

AP+ = {dP^/dk)Ak . (8) 

Furthermore, the period of a near-infrared field is about 
2.7 fs, over which the expectation value of the nuclear 
wavepacket changes by more than 2 a.u. (see. Fig. [3|. 
As long as the internuclear distance is not (indirectly) 
observed in an experiment, any interference effects will 
therefore smear out for near-infrared probe pulses with 
a few cycles of duration. Finally, strong-field ionization 
effects, such as rescattering ^ and multiple ionization 
bursts per half cycle of the field [43] , are other obstacles 
for the occurence of the double-slit interference effects in 
the partial differential ionization yields. 

Assuming an ideal scenario, in which the dissociating 
hydrogen molecular ion is ionized by an 800 nm laser 
pulse with just one cycle (total pulse length of 2.7 fs), we 
have performed numerical simulations for a pulse with 
a peak of the vector potential of Ao,i = 3 a.u. (about 
10^^ W/cm^ peak intensity). As shown in Fig. [t} we ob- 
serve a characteristic oscillation in the partial differential 
ionization yield at a specific value of the photoelectron 
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FIG. 7: (Color online) Results of numerical simulations for the 
partial differentail ionization yield of interaction with a 
1-cycle pulse at 800 nm. The photoelectron linear momentum 
was selected as ko = —1.36 a.u. and a detector resolution of 
Ak= 0.1 a.u. was considered. 



momentum of ko = —1.36 a.u.. The negative sign indi- 
cates that the electron is emitted in direction opposite to 
the laser polarization direction. We note that the yield 
for electrons emitted with the same absolute momentum 
value in the laser polarization direction, as well as the 
yields at different electron momenta, do not show any 
characteristic oscillatory behaviour. Thus, despite the 
fact that the oscillation in the specific yield in Fig. [T] 
appears to resemble the interference patterns, we expect 
that in general interference effects cannot be observed 
using a intense probe laser pulse at near-infrared wave- 
lengths. 



B. Imaging of nuclear wavepacket dynamics 

Now, we will analyze which information about the dis- 
sociating nuclear wave packet (velocity, spreading, etc.) 
can be retrieved from the double-slit interference pat- 
terns. For large internuclear distances we can well ap- 
proximate 

kR{r) + <!> ^ 2kvr kRc ^ ^ (9) 

where v is the velocity of the dissociating protons and 
kRc is a phase that depends on the initial internuclear 
distance and the initial acceleration due to the repulsive 
potential at small internuclear distances. Thus, the fre- 
quencies of the oscillations in Fig. |6] are approximately 
equal to 2kv and we can retrieve the velocity (or momen- 
tum) of the protons. From the results of our simulations 
for the 45 nm ionizing pulse, we read a frequency of 0.020 
a.u. (±0.001 a.u.), while the photoelectron momentum 
is k = 0.83 a.u. (±0.02 a.u.), obtained from the Fourier 
transform of the ionized photoelectron wave functions in 
the simulations. This gives rise to a velocity of v = 0.012 
a.u. (±0.001 a.u.). For the 22 nm ionizing pulse we 
obtain a frequency of 0.038 a.u. (±0.005 a.u.) and a 
momentum of k = 1.64 a.u. (±0.02 a.u.), resulting in a 
velocity of v = 0.011 a.u. (±0.001 a.u.). Both results are 
in good agreement with the velocity of v = 0.0104 a.u. 
(±0.0003 a.u.) obtained from the results for the expecta- 
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tion value of the internuclear distance (R), shown in Fig. 
ib). 

As outlined before, the decay of the oscihations in Fig. 
[6] depends on both the spread of the nuclear wavepacket 
in the internuclear distance, which increases with r, 
and the width of photoelectron momentum distributions, 
which remains constant with r in the present study. To 
estimate the decay we revise the formula for the double- 
slit interferences, given in Eq. ([7|), to account for the 
spreading as: 

/CO POO 
/ dkdR (10) 
-oo J — oo 

[1 + cos{kR + ^)] G{R - Ro, AR) G{k - ko, Ak) 



where G{a,b) is a Gaussian profile defined as 



G(a,6) = 



1 



27Tb 



e 



(11) 



with Rq = (R) being the mean value of the internuclear 
distance (see Fig. [2|b) for its dependence on r), AR is 
the width of the nuclear wavepacket (which depends on 
r as well), ko and Ak are the center and the width of 
the photoelectron momentum distribution, respectively 
(in an actual experiment, Ak may also be related to the 
resolution of the detector). Assuming that $ does not 
depend on the internuclear distance and the central pho- 
toelectron momentum is constant, we perform the inte- 
grations in Eq. (10) to get: 



(X 



1 ± — cos 
A 



kpRp 
A2 



2A2 



(12) 



where A = Vl + Ak'^AR^. kp and Ak are parameters, 
which can be retrieved in an actual experiment (here, 
we take these parameters from the results of our numer- 
ical simulations). Thus, we are left with two variables in 
Eq. (12), namely the internuclear distance Rq and the 



spreading of the nuclear wavepacket AR. As explained 
before, at large internuclear distances the Coulomb re- 
pulsion between the two protons becomes negligible and 
the protons propagate with constant momentum. In this 
regime we get 



Ro{t) ^ 2vr^Rc 



(13) 



and 



AR{r) ^ ^1 + 4r^Ak^/ml/{2Akp) + AR,, (14) 

where Akp is the spreading of the relative momentum of 
the protons and nip = 918 a.u. is the reduced mass of the 
two protons. AR, and R, are two terms which account 
for the initial nuclear spread, the initial internuclear dis- 
tance and the initial acceleration of the protons, where 
the linear temporal approximation fails due to the strong 
repulsive Coulomb potential. As shown before, v can be 
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FIG. 8: (Color online) (a) Ionization yields: Comparison 
between the results of numerical simulations (blue points) 
and the predictions of the analytical model, Eq. ( p^ , with 
(red dashed line) and without (green solid line) the cos-term, 
choosing v = 0.0111 a.u., k = 1.64 a.u., Ak = 0.084 a.u., 
Akp = 1.48 a.u., Rc = -0.45 a.u., ARc = -0.26 a.u. and 
$ = 1 a.u.. (b) Mean value of the internuclear distance: Com- 
parison between the results of numerical simulations (blue 
solid line) and the fit formula, Eq. (13) (red dashed line), (c) 
Spreading of the nuclear wavepacket: Comparison between 
the results of numerical simulations (points with error bars 
due to the grid resolution) and the fit formula, Eq. ( 14 ) (red 
dashed line). 



retrieved with high accuracy from the interference pat- 
tern. Assuming that Akp can be either observed in the 
experiment or approximately retrieved from the width 
of the photoelectron momentum distribution, we need to 
fit the two parameters R, and AR, to obtain (the mean 
value of) the internuclear distance and the width of the 
nuclear wave packet at the moment of ionization from 
the decay of the oscillations. 

In order to check our analytical model, we performed 
numerical simulations for the 22 nm ionizing pulse with 
a high resolution of the time delay {At = 0.04, = 15 
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as) (cf. points in Fig. [sj^a)). Choosing Rc = —0.45 a.u. 
and Ai?c = —0.26 a.u. we obtain a very good fit of 
the decay of the amplitudes of the oscillations (see green 
solid curves in Fig. [sj^a)). Using these parameters and 
choosing ^=1 rad, we were also able to fit the detailed 
oscillations very well (red dashed curve). As explained 
above, the fitting procedure allows us to further retrieve 
the mean value of the internuclear distance and the width 
of the nuclear wave packets. In Figs, [sj^b) and (c), we 
compare the results of these fits (red dashed lines) with 
the results retrieved from our numerical simulations (blue 
solid lines and points). The agreement is very satisfac- 
tory. Therefore, we can expect to be able to extract dy- 
namical information about the nuclear wavepacket from 
such interference patterns in an experiment. 

V. SUMMARY 

We have studied double-slit like interferences in the 
photoionization yields of a dissociating hydrogen molec- 
ular ion via numerical simulations of the time-dependent 
Schrodinger equation. To this end, we analyzed the dif- 
ferent steps in a pump-probe scenario, in which H2" is ex- 
cited by a pump pulse from the ground state to the first 
excited dissociative state and is then probed by a second 
time-delayed pulse, which ionizes the electron from the 
molecular ion as the ion dissociates. The results of our 
numerical simulations for a H2" model ion in its vibra- 



tional ground state show that the frequency of the pump 
pulse should be in resonance with the (jg — Gu transi- 
tion and should not be longer than about 5 fs in dura- 
tion (total pulse duration) in order to generate a well 
confined nuclear wavepacket, which is a requirement to 
observe interference effects in the dissociating molecular 
ion (as long as the internuclear distance of the protons 
at the moment of ionization is not observed). We have 
shown that the yields of ionization due to the interaction 
with a VUV probe pulse exhibit the expected double- 
slit interferences as a funtion of time delay between the 
two pulses. The amplitudes of the oscillation are found 
to decrease with increase of the time delay due to the 
spreading of the wave packet. By developing an analyt- 
ical model we were able to retrieve information about 
the nuclear wavepacket, namely its velocity, the mean 
value of the internuclear distance and the width of the 
wavepacket, from the oscillations in the ionization yields. 
We have tested the predictions of the analytical model 
against the results of numerical simulations and found a 
very satisfactory agreement. 
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